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Abstract 

We examine the B — )• D* form factor at zero recoil using a continuum QCD approach rooted 
in the heavy quark sum rules framework. A refined evaluation of the radiative corrections 
as well as the most recent estimates of higher order power terms together with more careful 
continuum calculation are included. An upper bound on the form factor of J^(l)^0.93 
is derived, based on just the positivity of inelastic contributions. A model-independent 
estimate of the inelastic contributions shows they are quite significant, lowering the form 
factor by about 6% or more. This results in an unbiased estimate 7^(1) ~ 0.86 with about 
three percent uncertainty in the central value. 



* On leave of absence from Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg 188300, Russia 



1 Introduction 



The determination of the CKM matrix element Vcb from exclusive decays has to rely on calcu- 
lations of the relevant form factors, which are usually defined as 

^(B^Diue) = —^\V,t>\'iMB + MnfMUco'-lf'ig{u)f (1.1) 

with u = V ■ v' = Ej-)(t) /Mp)(t) (in the B rest frame), and where -P(w) is a known phase space 
factor. Based on the normalization of the form factors at u = 1, Vcb is extracted from an 
extrapolation of the data to the non-recoil point. 

In the heavy quark limit, the normalization of the form factors J^(l) = ^(1) = 1 is given 
by heavy quark symmetry, and the main issue in the Vcb determination becomes a reliable 
calculation of the deviation from the heavy quark limit. The published extractions of Vcb along 
this route rely solely on the lattice calculations currently cited as [T] 

J^(l) =0.921 ±0.024 

g{l) = 1.074 ± 0.018 ± 0.016 (1.2) 

However, based on the dynamic heavy quark expansion in Minkowski space, it has been argued 
that larger deviations from the symmetry limit for J^(l) are natural in continuum QCD [2]; 
these have been further supported by the arguments [3J exploiting the relatively small kinetic 
expectation value fi^ extracted from the fits to inclusive B decays. The same line of reasoning 
led to a rather precise estimate of ^(1) in B^D [1], showing significantly smaller deviations 
from unity compared to Eq. (11. 2p . 

Recent years were of primary importance for heavy flavor physics. Along with advances 
in theory, many nontrivial nonperturbative predictions were verified with high precision, and 
heavy quark parameters experimentally extracted in accord with prior theoretical expectations; 
certain predictions were indirectly confirmed in dedicated lattice calculations. All this raised 
the credibility of OPE-based methods and favored an early onset of the short-distance expan- 
sion instrumental for high-precision predictions; confidence rose in the assumptions underlying 
dynamic treatment of the nonperturbative physics in heavy quarks. This progress warrants a 
critical re-examination of the form factors. The goal is to incorporate the accumulated knowl- 
edge and to shift the focus from merely establishing the scale of the deviations from the heavy 
quark symmetry limit towards obtaining a refined estimate with a motivated error assessment. 

In the present note we discuss the B ^ D* transition at zero recoil using a dynamic QCD 
approach inspired by the original treatment of the zero-recoil sum rules for heavy flavor tran- 
sitions. The details of the analysis will be presented in the extended publication [5]. 



2 Zero Recoil Sum Rule in QCD 

We consider the zero-recoil (^=0) forward scattering amplitude T^^{e) of the flavor-changing 
axial current cj'j^b off a i? meson at rest: 

T-(£) = j^'^j da^o e-^^«(*^^-*^°*^^)^(5|| tTc^,^,b{x) 6705c(O)|5) , (2.1) 
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Figure 1: The analytic structure of T^'^le) and the integration contour yielding the sum rule. 
Distant cuts are shown along with the physical cut. The radius of the circle is Em- 



where e is the excitation energy above Mo* in the B ^ Xc transition (the point e = corre- 
sponds to the elastic B — )■ D* transition). The amplitude T'''^{e) is an analytic function of e 
and has a physical decay cut at £:>0, and other distant singularities. The analytic structure of 
T^^{e) is shown in Fig. [H 
The contour integral 

/o(^A/) = -7^ lT^'{e)de (2.2) 
Zm J 

|e|=eAf 

with the contour running counterclockwise from the upper side of the cut, see Fig. [H leads to 
the sum rule involving -7^^(1). Using the analytic properties of T^'^{£) the integration contour 
can be shrunk onto the decay cut; the discontinuity there is related to the weak transition 
amplitude squared of the axial current into the final charm state with mass Mx = M£,*+e. 
Separating out explicitly the elastic transition contribution B^D* at e = we have 

/o(£m) = J^^(l) +Winel(£M), Wi„cl(£Af) = ^ j dlSC T^' {e) ds , (2.3) 

e>0 

where Wmei(^A/) is related to the sum of the differential decay probabilities into the excited 
states with mass up to Mr)*+£Ai in the zero recoil kinematics. 

The OPE allows us to calculate the amplitude in (12.11) - and hence Io{sm) - in the short- 
distance expansion provided \e\ is sufficiently large compared to the ordinary hadronic mass 
scale. It should be noted that strong interaction corrections are driven not only by \e\, but also 
by the proximity to distant singularities. Therefore, Em cannot be taken too large either, and 
the hierarchy e^f <^2mc has to be observed. 

The sum rule Eq. (12. 3p can be cast in the form 

J^{1) = VH£M)-WineiieM) (2.4) 

which is the master identity for the considerations to follow. Since Wmei(£^Af) is strictly positive, 
we get an upper bound on the form factor 

^(1) < Vloi^M) (2.5) 

which relies only on the OPE calculation of Jq. Note that this bound depends on the parameter 
Em, while (12.41) is independent of Em since the dependence in /q and Winei cancel. Furthermore, 
including an estimate of Winei(£^Af) we obtain an evaluation of J^(l). 

The correlator in (12. ip can be computed using the OPE, resulting in an expansion of T^'^Ie) 
in inverse powers of the masses rric and nib- This results in the corresponding expansion of 
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IqIem)- This OPE takes the following general form 

loieAi) = C/' (^M, /^) + V Ckieu, /i) — - — 3 (2.6) 

k "^Q 

= ^A^^^i^M^f^) — ^l/m^(^M, At) — ^l/m|(^M,/^) " (^^M, /^) " ••• 

where Ok are local 6-quark operators b...b of increasing dimension dk > 5, Ck{lj) are Wilson 
coefficients for power-suppressed terms, and ^^^^^^ is the short- distance renormalization (corre- 
sponding to the Wilson coefficient of the unit operator), which is unity at tree level. We have 
also introduced a Wilsonian cutoff /x used to separate long and short distances. The complete 
result, of course, does not depend on n since the /i-dependence cancels between the Wilson 
coefficients and the matrix elements of the operators. At tree level A does not depend on sm- 
The choice of /i is subject to the same general constraints as that of Em-, and it is therefore 
convenient to choose ^ = eM- 



2.1 Perturbative corrections 

The perturbative renormalization ^^^^{fJ-) can be expanded in power series in as- We use the 
Wilsonian OPE and benefit from well-behaved perturbative series for ^"^^^{n)- The exact form 
of the perturbative coefficients depends on the definition chosen for higher- dimension operators; 
we adopt the often used kinetic scheme [B]. 

In one-loop perturbative calculations there is a simple connection between the normalization 
point of the heavy quark operators in the kinetic scheme and the hard cut-off on the gluon 
momentum in the diagram. This allows to obtain the analytic expression for ^^'^'^*(/i) to this 
order even without explicit calculation of the Wilson coefficients Ck in Eq. (12. 6p . The expression 
is rather lengthy and will be presented in Ref. |S]. By the same trick one also obtains all higher- 
order BLM corrections by performing the one- loop calculations with massive gluon [7j. 

A similar argument does not apply to non-BLM corrections starting where £M-dependence 
of ^"^^^ has to be determined expanding in l/m^; for C(a^) corrections this was done in Ref. ([8]) 
through order I/ttiq. The corresponding coefficient was found to be small numerically, which 
suggests that omitted terms oc a^ef^/mg and higher should not produce a significant change. 

Perturbative corrections to ^"^^^{em) appear to be small for practical values of Em between 
0.6GeVand 1 GeV. Taking, for instance, £:M = 0.75GeV, mc = 1.2 GeV, mb = 4.6GeV, as{mb) = 
0.22 we get the numeric estimates at different orders 

= 1 - 0.019 + (0.007 - 0.004) + 0.0045 + ... (2.7) 

Here the first term is the tree value, second is 0{as) evaluated with as = 0.3, the next pair of 
values show the shift upon passing to the C(a^) order (positive for the BLM part and negative 
from the non-BLM contribution); the last term shows P^a^ term as an estimate of even higher- 
order perturbative corrections. Fig. [2] shows the dependence on Em of these predictions for 

^Jc^^^■ In particular, taking the full two-loop result as the central estimate we find 

y^^^''*(0.75 GeV) = 0.985 ± 0.01 ; (2.8) 

we will use Em = 0.75GeV in the subsequent discussion. We emphasize that the numeric stability 
applies only to the perturbative renormalization factor in the Wilsonian OPE, and the quoted 
values refer to the specific renormalization scheme (kinetic) adopted in the analysis. 
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Figure 2: Left: 



(^^^^ to order Og (blue), including /3oa'f (green), full 0{a'j) (red) and including 

PqU^ (magenta), assuming a^^(mb) = 0.22, mc = 1.2GeV and mfe = 4.6GeV. 
Right: Upper bound (j2.5p on depending on em, with or without PqU^ term. A 

fixed A = 0.11 is used; assuming the perturbative evolution of A with em would 
flatten the dependence. 



2.2 Power corrections 

The leading power corrections to Iq were calculated in Refs. |9l [10] to order l/rrig and to order 
l/rriQ in Ref. [H] and read 



^l/m2 



f^G _|_ f^G 



3m^ 



+ 



Arrfi 



+ 



1 



+ 



12mh 



mcTTii, mi 



+ Z3 (PD + Pis) 



(2.9) 



The nonperturbative parameters /x^, /i^, p|, and p^^. all depend on the hard Wilsonian cutoff. 
The ^M-dependence of ^^^^^^ is linked to the power-like scale dependence of the nonperturbative 
matrix elements through their mixing with lower-dimension operators. 

In the numerics we use the values /i^(0.75 GeV) =0.4 GeV^ pl)(0.75 GeV) =0.15 GeV^ while 
for the quark masses mc = 1.2 GeV and m;, = 4.6 GeV (the scale dependence of the latter plays 
a role here only at the level formally beyond the accuracy of the calculation). The dependence 
on p^ and on p\g is minimal and their precise values do not matter; we use for them 0.3 GeV^ 
and — 0.12GeV^, respectively. We then get 



Ai 



0.091, 



0.028 



(2.10) 



If we employ the values of the OPE parameters extracted from a fit to inclusive semileptonic 
and radiative decay distributions [I2l |13] , we find a consistent result 



Ai/„2 + Ai/„3 = 0.11 ± 0.03 , 



(2.11) 



An important question is how well the power expansion for the sum rule converges. Recently, 
the OPE for the semileptonic S-meson structure functions has been extended to order I/^q 
and l/mg [T5] . Combined with the estimates (THj of the corresponding expectation values 
discussed in Ref. |15|, this leads to 



Ai/^4 ^ -0.023, 



A , r ~ . 



-0.013. 



(2.12) 



We then observe that the power series for Iq appears well-behaved at the required level of 
precision. For what concerns the loop corrections to A, the 0{as) correction to the Wilson 
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coefficient for the kinetic operator in Eq. fl2.9p was calculated in Ref. ^ and turned out nu- 
merically insig nificant. At C(l/m^), even if radiative corrections change the coefficient for the 
Darwin term by 30% the effect on the sum rule would still be small. 

Taking into account all the available information, our estimate for the total power correction 
at 5M = 0.75GeV is 

A = 0.105 (2.13) 

with a 0.015 uncertainty due to higher orders. On theoretical grounds, larger values of /i^ and/or 
are actually favored; they tend to increase A. Combining the above with the perturbative 
corrections we arrive at an estimate for Iq and, according to Eq. ( 12. 5p at a bound for the form 
factor, which in terms of central values at eA./ = 0.75 GeV is 

J^(l)<0.93. (2.14) 

As stated above, the upper bound in Eq. (12.51) depends on Em, see Fig. [21 becoming stronger 
for smaller Em- It is advantageous to choose the minimal value of Em for which the OPE-based 
short- distance expansion of the integral (12. 2p for Io{em) sets in. This directly depends on 
how low one can push the renormalization scale /i while still observing the expectation values 
actual /i-dependence in the kinetic scheme approximated by the perturbative one. Since in this 
scheme /i^(/i) >/iG'(/^) holds for arbitrary /x, in essence this boils down to the question at which 
scale /imin the spin sum rule and the one for /i^, get approximately saturated, e.g. /^^(/imin) — 
0.3 GeV^. The only vital assumption in the analysis is that the onset of the short-distance 
regime is not unexpectedly delayed in actual QCD and hence does not require Em > 1 GeV. 
This principal question can and should be verified on the lattice. This will complement already 
available evidence from preliminary lattice data [LTj as well as from the successful experimental 
conffimation [TS] in nonleptonic B decays of the predicted 3/ 2- dominance. 



2.3 Estimate of the inelastic contribution Winei 

On general grounds Winei is expected to be comparable to the power correction A considered 
above. To actually estimate it we consider another contour integral 

/i(£A/) = -^ / ^"(^)^d^ (2.15) 

|e|=eA/ 

for which we can write 

^inel(eM) = (2.16) 

where £ is the average excitation energy (it depends on em)- The integral is expected to be 
dominated by the lowest radial excitations of the ground state, with e;^erad~700MeV. 
-^i(^m) can also be calculated in the OPE flO]; the result including l/m^ terms reads 

^ -(pIg + Pa) ^ -2pL - pIg ^ pL + pIg + Ps + Pa M ^ 2 ^ J_ 
3m^ 3mcmb 4 V^c Snicmb ml 

where the non-local zero momentum transfer correlators p^^, p^Q, p| and p\ are defined in |10] . 
They can be estimated along the lines described in [15], based on saturation by the appropriate 
intermediate states. We shall defer the details of this estimate to [5]. Here we note that only 
the ffist term survives in the BPS limit where it is positive; the second and third terms are 
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of first and second order in the deviation from the BPS hmit, respectively. The last term is 
positive being the correlator of two identical operators h{(T'K)%. Since the middle term comes 
with a small coefficient l/3mcmfe, the expression has only a shallow minimum where the first 
term is decreased by less than 10%; the sum of the last two terms becomes larger than that 
only if it is positive. On the other hand, the combination 

- {pIg + Pa) + pIs (2.18) 

determines the hyperfine splitting to order I/ttiq and can be constrained from the observed 
masses of 5^*^ and D^*^ mesons. We finally obtain 

, , , . 0.48 GeV^ 

with some uncertainty from perturbative corrections, implying for e = Crad — 700 MeV 

Winel ~ — > 0.13 . (2.20) 

Crad 

This estimate is derived at leading order in l/mq and may be corrected by higher-order 
terms by as much as 30%. We observe that i^inei is similar in size to the power term in /q and 
even exceeds it. Using (12.201) at face value we arrive at our estimate for the expected value of 
the form factor 

J^(l) < 0.86. (2.21) 

The quasi- resonant states are expected to dominate ifinei(yu) at intermediate yU ~ 1 GeV; 
the continuum contribution to w\-ae\{p) is parametrically l/A^^c-suppressed and usually smaller. 
The D^*^7r continuum can independently be evaluated in the soft-pion approximation. It turns 
out that numerically the dominant effect originates from the heavy quark symmetry breaking 
difference between the B*Btt, D*Dtt and Z}*Z}*7r-couplings which until recently has not been 
accounted for in this context jl| although is expected to be significant [19]. The result is shown in 
Fig. [3l it depends on the upper cutoff in the pion momentum marking the end of the soft 
continuum domain for D^*^7i, presumably somewhat below erad- We expect about 5% combined 
yield for Tjj*+ = 96 keV, i.e. about a third of the overall Wi^ei in E<1- (I2.20p in accord with the 
1/Nc arguments. This contribution alone would lower the upper bound in Eq. (I2.14p by 0.025. 



3 Conclusions 

The direct OPE-based l/rriQ expansion of the zero recoil form factor which we have analyzed 
through the sum rule for the correlator of the zero-recoil axial currents yields the unbiased 
estimate J^(l) ~ 0.86, and suggests a lower bound J-'{1) < 0.93. All the values refer to pure 
QCD form factors and exclude possible electromagnetic effects, in particular the universal 
short- distance semileptonic enhancement factor of 1.007. 

Since a charm mass expansion is involved, the precision of our estimates is limited. Accuracy- 
wise two aspects should be kept distinct. On the one hand, the heavy quark parameters enter 
our analysis, and in particular /x^, and m^. Their determination from inclusive measure- 
ments will soon be improved thanks to refined theoretical calculations. On the other hand, 

^The D*D*Tr channel has not been previously considered while generally required by the heavy quark sym- 
metry. 
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Figure 3: Non-resonant D-k and D*tt contribution to winci at (7D*£i7r = 4.9 GeV corresponding to 
rr,*+=96keV. The plots show, from bottom to top, Sbibj^^i q.S, 0.6 and 0.4, and 1bibj^ = i q.S 
and 0.6, respectively 



even if these QCD parameters were accurately known, sharpening the upper bound and the 
estimate of would require additional nonperturbative input. To some extent the uncer- 
tainties may be reduced by dedicated measurements in semileptonic decays or by direct lattice 
calculation of the relevant heavy quark parameters. In our opinion, the latter can be done in a 
straightforward way, which will be discussed elsewhere. 

In our numerical estimates we used the lowest values of /x^ and consistent with the 
theoretical bounds; larger values would typically further lower both the upper bound and the 
central expectation for J-'(l). We estimate the perturbative and non-perturbative uncertainties 
in our approach to be about iO.Olpt and ±0.02np (iO.Olnp for the upper bound), mostly related 
to the unavoidable truncation of the l/nic expansions. The uncertainty in the prediction is 
dominated by winei, is not symmetric, and allows for larger deviations towards lower values of 
^(1). 

There is an alternative way of determining \Vcb\i using inclusive semileptonic decays. It 
rests on the heavy quark expansion in l/m^, and on the OPE which, in this case, should be 
considered a mature tool. Recent estimates of higher orders in l/mq as well as in as did not 
show signs of failure of the heavy quark expansion in this application. Hence one can employ 
the value of Vcb = (41.54 ± 0.44 ± 0.58) x 10^^ from the inclusive decays [12] and determine the 
form factors from the data on B D^*^^ decays; this yields 

J^(l) = 0.86 ±0.02 (3.1) 
e;(l) = 1.02 ±0.04 (3.2) 

The value for ^(1) is in a perfect agreement with the HQE prediction [3|; J^(l) hits the cen- 
tral value of our estimate. Both form factors in Eq. (13. ip are noticeably below the quoted 
lattice values Eqs. (11.21) : this explains the observed tension between the inclusive and exclusive 
determinations of |Vc6|. 

Our analysis suggests that the lattice determinations of both form factors are systematically 
high. The central values of the currently quoted lattice J^(l) is very close to our upper bound. 
Lattice calculations - provided they accommodate the experimentally measured 5-meson ex- 
pectation values - seem to imply an extremely small inelastic contribution, which is a priori 
unnatural. Moreover, it appears in contradiction with the large non-local correlators encoun- 
tered to order l/wig. The estimate of the non-resonant soft-pion D^*^n rates also confirms 
this assessment. We conclude that values for J-'(l) in excess of 0.9 would be consistent with 
unitarity and the short- distance expansion of QCD only with rather contrived assumptions. 
Values of J^(l) larger than 0.93 should be viewed as violation of unitarity assuming that usual 
short- distance expansion in QCD works. 
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With the unbiased estimate based on the lower-end /i^ and p|, yielding J-'(l) ~ 0.86 - or 
even smaller for larger /i^, pf, - we find a surprisingly good agreement (probably, somewhat 
accidental) between the exclusive and inclusive approaches. In view of the inherent l/rric 
expansion for and of the approximations used for proliferating hadronic expectation values, 
matching the theoretical precision already attained for Vcb from the inclusive fits does not look 
probable. Nevertheless, additional experimental and/or lattice input would make the suggested 
3% uncertainty interval more robust. 

We estimate that the contribution of excited radial states with mass below 3 GeV constitutes 
about 10% or more of the total yield. This refers only to the zero-recoil kinematics and only 
to the axial current; it does not include P-wave states. This fraction may even be larger when 
applied to the full phase space. This suggests that the observed 'broad state' yield in this mass 
range routinely attributed to the | P-wave excitations is actually dominated by states with 
different quantum numbers, thus resolving the '| > | puzzle'. The suppression of the broad 
P-wave yield was predicted based on the spin sum rules and confirmed indirectly in nonleptonic 
B decays [18]; a recent discussion can be found in Ref. [20] . 
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